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Lautzenheiser  and  Edward  Sable,  as  well  as  the  National  Climatic  Center's  spec- 
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student  assistants  under  Sister  Leonards  Burke,  Regis  College,  Weston,  MA.  The 
third  example  was  initially  undertaken,  frankly,  with  doubts  as  to  its  feasibility. 
But  the  idea  was  conceived  and  encouraged  by  my  associates,  Messrs.  Grantham 
and  Lund,  who  also  undertook  much  of  the  preliminary  planning.  Mr.  Grantham 
instructed  the  National  Climatic  Center  assistants  in  Asheville  on  the  data  collec- 
tion. Funds  for  the  radar  study  were  provided  by  the  United  States  Air  Force. 

The  function  of  the  AFCRL  Computation  Branch's  Analysis  and  Simulation  Unit 
made  this  study  possible  through  the  excellent  programming  by  the  contractor, 
Boston  College,  on  the  AFCRL  computer  installation. 
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Areal  Coverage  Estimates 
by  Stochastic  Modelling 


1.  INTRODUCTION 

At  a point  of  observation,  how  representative  is  the  weather  with  respect  to 
the  surrounding  area?  With  the  passage  of  a storm,  what  fraction  of  a city  or 
state  will  have  received  over  25  mm  of  rainfall,  a trace  or  more,  or  will  have 
been  missed  by  the  storm?  What  is  the  probability  that  a satellite  camera  will 
record  a clear  view  of  a given  area  of  the  globe,  or  a percentage  of  the  area? 
These  and  many  more  questions  on  the  horizontal  persistence  or  variability  in  the 
weather  are  similar  to  the  questions  asked  with  regard  to  temporal  events.  But 
whereas  the  latter  have  been  the  subject  of  many  successful  studies  and  papers, 
the  problem  of  areal  coverage  has  remained  challenging. 

Previous  studies  of  the  areal  extent  of  the  weather,  for  the  most  part,  have 
yielded  essentially  spatial  correlation  coefficients,  but  there  have  been  notable 
exceptions..  Riedel  et  al*  worked  with  integrated  rainfall  amounts  in  areas  of  given 
size,  making  estimates  of  the  maximum  possible  precipitation  falling  in  periods  of 
6 to  48  hours.  They  obtained  the  variation  of  extreme  total  rainfall  for  areas 


(Received  for  publication  6 July  1976) 

1.  Riedel,  J.F. , Appleby,  J.F. , and  Schloemer,  R.  W.  (1956)  Seasonal  Variation 
of  the  Probable  Maximum  Precipitation  East  of  the  105th  Meridian  for  Areas 
from  10  to  1000  sq,  miles  and  Durations  of  6,  12,  24,  and  48  Hours.  U.  S. 
Dept,  of  Commerce,  Hydrometeorological  Report  No.  33,  Washington,  D.  C. 
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ranging  from  26  to  2600  km  . Among  studies  of  storms.  Court  reviewed  and 


tested  models  of  rainstorms  that  were  assumed  to  be  symmetrical  around  a storm 
center.  Likewise  the  assumption  of  symmetry  was  basic  to  the  models  of  Briggs' 
and  Roberts.  The  latter  defined  a scale  measure,  on  the  order  of  80  to  320  km, 
for  the  radial  extent  of  precipitation,  which  he  derived  from  the  frequency  of  the 
event  at  a single  station  and  its  probability  of  occurrence  somewhere  in  a dense 


network  of  stations. 

More  recently,  Sims  and  Jones  surveyed  the  rainfall  rates  averaged  over 
lines  10  to  62  km  long  in  two  special  field  projects,  one  in  central  Illinois,  the 
other  in  Florida.  The  latter  project  was  part  of  the  well-known  1948  Florida 
Thunderstorm  Project.  Sims  and  Jones'  recourse  to  the  data,  now  more  than 
25  years  old,  is  indicative  of  the  infrequency  of  experimental  surveys  of  areal 
coverage  on  a mesoscale;  however,  if  experimental  efforts  have  been  few,  models 
of  areal  coverage  are  virtually  nonexistent. 

This  paper  makes  no  effort  to  describe,  or  map,  actually  occurring  details  of 
weather  phenomena;  rather,  its  purpose  is  to  estimate  the  probability  of  events 
occurring  along  lines  of  travel  or  within  specified  areas  or  fractions  thereof.  In 
contrast  with  Riedel's  purpose  to  present  integrated  rainfall  amounts  in  a given 
area,  the  effort  of  this  paper  is  centered  on  the  maximum  or  minimum  value  of  a 
meteorological  element,  either  in  a given  area  or  along  a line  of  travel.  The 
archived  data  are  used,  not  as  comprehensive  climatology  in  themselves,  but  as  a 
guide  to  modelling.  Additionally,  a model  is  developed  to  yield  answers  on  the 
fractional  cover  of  an  area. 

2.  APPRO \CH 

The  problem,  at  the  outset,  is  to  devise  a procedure  that  will  produce  a sim- 
ulated horizontal  synoptic  field  stochastically,  without  recourse  to  physical  laws  or 
dynamics.  Once  this  is  accomplished,  the  procedure  can  be  repeated  many  times 
to  produce  a large  set  of  synoptic  situations,  thus  comprising  a sample  that  can  be 
surveyed  for  spatial  correlation  or  association  of  events. 


2.  Court,  A.  (1961)  Area-point  rainfall  formulas,  J.  Geophys.  Res. , 66:1823- 

1831. 

3.  Briggs,  J.  (1972)  Probability  of  aircraft  encounters  with  heavy  rain, 

Meteorol.  Mag. , 101:8-13. 

4.  Roberts,  C.  F.  (1971)  A note  of  the  derivation  of  a scale  measure  for  precipi- 

tation events.  Monthly  Weather  Review,  99:873-876. 

5.  Sims,  A.  L.  and  Jones,  D.  M.A.  (1975)  Frequencies  of  short-period  rainfall 

rates  along  lines,  J.  Appl.  Meteorol. , 14:170-174. 
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It  is  assumed  that  at  any  single  point  in  the  atmosphere  the  process  of  change 
of  a weather  element  is  stationary  stochastic;  that  is,  it  has  a partially  random 
pattern  of  change,  and  the  pattern  does  not  change  with  time.  Further,  it  is  sup- 
posed that  any  meteorological  quantity  <x)  can  be  transformed  into  the  normalized 
variable  (y)  in  a one-to-one  correspondence.  Conventionally  the  normalized  vari- 
able is  written  as  y(0,  1)  meaning  that  it  has  a mean  of  0.  0 and  a variance  of  1.0 
with  the  bell-shaped  distribution  implied.  The  transformation  (x  - y)  might  be 
accomplished  through  algebraic  or  analytical  expressions,  but  not  necessarily. 
Frequently,  it  is  sufficient  to  obtain  the  distribution  of  y graphically,  or  tabularly, 
through  the  cumulative  probability; 

y 

P(x>  P(y)  1/  j exp  (*?2/2)  d?  , d) 


where  the  symb  P is  used  for  cumulative  probability  or  the  probability  that  the 
value  in  parenthesis  is  not  exceeded. 

The  transformation  (x  - y)  is  possible  for  weather  elements  whether  they  are 
continuous,  like  air  temperature,  or  classified  in  categories,  like  sky  cover  (clear, 
scattered,  broken,  or  overcast).  In  one  example  (Figure  1)  the  transformation  of 
the  Minneapolis  January  hourly  temperature  (T)  into  y(0,  1)  is  straightforward. 

The  ordinate  scale  is  the  cumulative  probability  P(T)  5 P(y),  and  the  abscissa  is 
the  temperature  (T).  The  value  of  y(0,  1)  is  shown  alongside  P(y),  but  is  available 
in  almost  any  textbook  in  statistics.  Through  the  plotted  curve  (Figure  1)  the 
graphical  transformation  of  T into  y(0,  1)  is  accomplished. 

In  another  example  (Table  1)  the  sky  cover  in  January  at  midnight  in 
Minneapolis,  Minnesota,  has  been  reported  in  frequencies  of  33  percent  for  clear, 
12  percent  for  scattered  clouds,  11  percent  for  broken  clouds,  and  44  percent  for 
overcast.  From  a table  of  the  normal  distribution,  the  variable  (y)  takes  on  values: 
«-0.43  for  clear  sky,  between  -0.43  and  -0.  13  for  scattered  clouds,  between 
-0.  13  and  0.  16  for  broken  clouds,  and  >0.  16  for  overcast. 

At  each  station  at  each  time  of  day  in  each  month  or  season,  an  element  like 
temperature  will  have  its  own  characteristic  distribution  of  values.  But  all  trans- 
formed variables  (y),  from  all  stations,  will  have  a common  distribution  with  a 
mean  of  0.0,  a variance  of  1.0,  and  the  Gaussian  density  function.  The  goal  of  this 
paper  is  to  devise  a stochastic  process  with  the  normalized  variable  generated 
everywhere.  Between  the  values  of  y at  any  two  stations  there  will  be  a correlation 
that  generally  will  decrease  with  increasing  distance  between  stations.  Within  a 
given  area  (circle,  square,  or  other  shape)  there  will  be  a set  of  values  that  will 
have  a minimum  value  and  a maximum,  also  a highest  minimum  (or  lowest 
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Figure  1.  Plot  on  Normal  Probability  Paper  to 
Illustrate  the  Transformation  of  a Variable  (T) 
into  the  Normalized  Variable  y(0,  1)  Through  the 
Cumulative  Frequency  of  Values 


Table  1.  Example  of  the  Transformation  of  Categories  of  a Weather 
Condition  Into  Intervals  of  the  Normalized  Variable  (y).  The  exam- 
ple is  for  sky  cover  at  Minneapolis,  MN. , midnight,  January 


Normalized 
Variable  (y) 


Broken 


maximum)  value  for  a specified  fraction  of  the  area.  Each  minimum  or  maximum 
can  be  transformed  back  to  the  original  weather  element  through  reverse  applica- 
tion of  Eq.  (1). 


Suppose,  initially,  random  normal  numbers  (y  ) are  assigned  to  all  points  of 
a square  grid  with  ^ fine  mesh  (Figure  2).  At  a given  point  (i,  j)  suppose  the  num 
ber  yQ(i,  j)  is  replaced  with  a weighted  sum  of  the  surrounding  random  numbers. 
If  the  mesh  is  infinitesimally  fine'  and  the  area  infinite,  then  the  new  number  is 


Figure  2.  A Grid  of  Points  to  Which  Random 
Numbers  Are  Initially  Assigned.  The  random 
value  yQ(i  + ?,  j + rj)  is  at  distance 

R = V (€2  + r/2)  from  (i,  j).  The  value  y(i,  j) 
is  a weighted  average  of  the  random  numbers 
surrounding  (i,  j)  within  radius  Rq.  The  cc 
between  y(i,  j)  and  y(i  + k,  j + f ) decays  with 

the  distance  s = V(k2_+  {2) 


where  <?,  n)  are  coordinates  of  <i  + S,  j + p)  relative  to  (i,  j).  R is  the  radial 
distance,  given  by 


i. 


R” 


2 


+ 


p 


2 


(3) 


f ( R ) is  the  weighting  function  that  generally  should  decrease  with  increasing  dis~ 
tance  (R),  and  D is  a constant,  to  be  determined  such  that  the  variance  of  y(i,  j) 
becomes  1.0. 

At  another  point  (i  + k,  j + I ),  at  a distance  » f.om  (i,  j)  given  by 


the  random  number  is  replaced  with 


oo  oo 

y(i  + k,  j + /)  = -g-  J J F(R')  • yQ(i  + k + f',  j + / + rj')  dC  dp'  , (5) 

”0C  “00 


where  (£',  p')  are  coordinates  relative  to  (i  + k,  j + f)  and 

R'2  = C'2  + p'2  . (6) 

The  correlation  coefficient,  hereinafter  shortened  to  cc  and  symbolized  asp  , 
between  y(i,  j)  and  y(i  + k,  j + t ),  is 

p(s)  = E { y( i,  j)  • y(i  + k,  j + f)}  , (7) 

From  (2),  (5),  (7)  it  is  seen  that  p(s)  can  be  given  in  terms  of  an  infinite  number 
of  the  products  of  pairs  of  initial  random  numbers.  But  the  expected  product  of 
any  two  random  numbers  is  zero,  unless  they  coincide  in  which  case  their  expected 
product  is  1.0.  Hence,  using  (3),  (4),  (6) 

00  00 

p(s)  = — J J F(R)  • F(R')  d?  dp  , (8) 

“00  “00 
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when 


' 


k 


S'  = S - k , (9) 

n'  = n - t 

For  3 = 0,  p (s)  = 1,  hence 
00  00 

D 2 = / / F2(R)  dS  dr?  (10) 

-oo  -oo 
or 

00 

D2  = 2tt  j R • F(R)  dR  . (11) 

o 

In  practice,  for  a grid  that  is  less  than  infinitesimally  fine,  (2)  is  approxi- 
mated by 


y(i,  j)  = -jj-  S F(R)  • y (i  + S,  j + n)  , (12) 

where  (S,  rj)  are  given  only  in  whole  units  and  their  vector  sum  R s Rq,  where  Rq 
is  large. 

For  a finite  grid  (3)  holds  for  R,  (6)  for  R',  and  (8)  and  (10)  are  approximated 

by 


p(s)  = £F(R)  • F(R')  , (13) 

D2 

where 

D2  = 2 F2(R)  . (14) 

n 

If,  initially,  random  numbers  (x)  with  a rectangular  distribution  between  0.  0 
and  1.0  are  assigned  to  the  grid  points,  their  mean  will  be  1/2  and  variance  equal 
to  1/12.  Equation  (12)  is  given  alternatively  by 
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L 


i 


In  a Monte  Carlo  exercise  this  form  of  the  derivation  of  y-values  from  a set  of 

random  numbers  (x  ) becomes  a valuable  tool, 
o 

So  far,  the  model  has  been  considered  in  general  terms.  To  proceed  further, 
the  form  of  the  weighting  function  F(R)  must  be  specified.  In  the  following  sub- 
sections several  possibilities  are  considered,  including  two  of  specific  interest. 

:i.  1 ft  hen  Weighting  Function  Decreases  FxpontialU  with  Distance 
Let  the  weighting  function  F(R)  be  given  by 

F(R)  = exp (-Rb/A)  , (16) 

where  b is  kept  within  the  limits 

0 < b £ 2 . (17) 

The  parameter  A strengthens  or  weakens  the  horizontal  persistence.  Equation  (11) 
becomes 

m!r(A/2)m  • r(m)  , (18) 

where  m = 2/b  and  F(m)  is  the  gamma  function  of  m.  For  the  cc  between  points 
separated  by  s,  (8)  becomes 

00  oo 

p(s)  - -Tj-  J J exp  jj-(Rb  + R'b)/Aj  d?  dp  . (19) 

^ -00  -00 

Except  when  b = 2,  attempts  at  rigorous  analysis  of  (19)  failed  to  produce  a 
simple  analytic  solution  for  p(s).  It  was  necessary  to  approximate  p(s)  by  numer- 
ical process  on  an  electronic  computer.  This  was  done,  using  the  Gaussian- 
Laguerre  24 -point  double-precision  formula^  with  values  of  b ranging  from  0.  25 
to  1.5.  The  expertise  was  provided  by  Dr.  Richard  Dougherty,  RDP  Corporation, 
Bedford,  Massachusetts,  serving  as  a consultant  under  USAF  contract.  His 

(L  Stroud,  A.  H.  and  Secrest,  D.  (1966)  Gaussian  Quadratic  Formulas.  Prentice- 
Hall]  Englewood  Cliffs,  New  Jersey. 
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results  showed  that  the  weighting  function  (16)  produces  a field  in  which  the  cc 
decreases  exponentially  as  the  square  of  the  distance  for  short  distances,  then 
decreases  more  slowly  with  increasingly  large  distances.  Such  a correlation  field 
does  not  conform  to  most  previous  concepts.  For  the  present  the  model  of  (16)  has 
been  set  aside,  except  when  b = 2,  as  in  the  following  section. 

3.2  Model  V When  Weightin':  Function  is  Made  Proportional  to  the  Bivariate  Normal 
Density  Function 

Let  the  weighting  function  F(R)  be  given  by 

F(R)  = exp (-R2/A)  , (20) 

where  A again  is  a parameter  that  strengthens  or  weakens  the  horizontal  persist- 
ence. Then  (11)  becomes 

D2  - trA/2  (2D 

and  (8)  resolves  and  becomes 

p( s)  = exp  (-s2/2A)  . (22) 

Thus,  the  cc  is  found  to  decay  exponentially  with  the  square  of  the  distance  between 
stations. 

It  is  convenient  to  set  A = 50,  making  the  cc  equal  very  nearly  to  0.  99  for  unit 
distance  (s  = 1). 

To  conform  to  previous  conceptions,  it  is  desirable  to  find  a field  in  which  the 
cc  decreases  exponentially  with  distance  (not  squared),  at  least  for  short  distances. 
Additionally  there  would  be  an  advantage  to  a model  in  which  the  cc  drops  to  zero 
eventually  in  a large,  but  finite,  distance.  Thip  is  accomplished  in  Model  B. 

3.3  Model  B:  When  the  Weight  is  Made  Uniform  Within  a Circular  Area  Surrounding  the  Station 

Let  the  weighing  function  become  simply 

F(R)  = 1 for  R < Rq  , (23) 

= 0 for  R > R , 

o 
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where  R is  hereafter  termed  the  Radius  of  Influence  (Figure  2).  This  concept  of 
the  weighting  function  may  be  resisted  at  first,  but  it  resembles  the  concept  of  the 
rain  cell  of  Briggs3  and  Rq  resembles  the  "scale  measure"  of  Roberts.  The  value 
of  y(i,  j)  at  point  or  C2  (Figure  3)  is  a simple  average  of  many  uniformly  spaced 
random  numbers  inside  the  circular  region  of  radius  RQ  that  surrounds  the  point, 
with  a suitable  value  for  D to  make  the  variance  of  y(i,  j)  equal  to  unity.  The  cc 
between  the  two  normalized  values  at  Cj  and  Cg  is  positive  only  if  there  is  an 
overlap  of  the  two  circles,  and  is  directly  proportional  to  the  area  of  intercept. 

For  the  distance  s between  and  C^,  OS),  (11),  and  (4)  yield 


= 0 for  ct  2 1 , 


where 

a = s/2R 


1 

1 c,' u 

! \ 

Lj! - 

' 

w 

^ ) 

Figure  3.  Illustration  of  the  Radius  of  Influence  (RQ). 
In  Model  B the  intersection  of  the  two  circles  is  a 
direct  measure  of  the  cc  between  values  of  a weather 
element  at  Cj  and  C2>  decreasing  with  increasing  dis- 
tance (s) 


For  a small  distance  s,  the  Taylor  expansion  of  terms  of  (24)  is  approximated  by 

p (3)  * exp  (-2s/irRo)  . (26) 

That  is,  the  cc  decreases  exponentially  with  short  distances  (not  squared)  between 
stations.  The  departure  of  (24)  from  (26)  is  less  than  0.  1 for  s < 0.  8Ro>  But  as 
mentioned  in  Section  3.  2,  the  decrease  of  cc  to  zero  at  distance  2Rq  should  be  an 
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advantage  considering  that  the  cc's  reported  by  most  authors  do  drop  to  zero  at 
some  measurable  distance  between  stations. 


It  is  convenient  to  set  R^  = 64  to  give  the  cc,  by  (24),  very  nearly  equal  to 


0.99  for  a unit  distance  (s  - 1). 


t.  SCALE  DISTANCE 


In  Section  3,  distances,  symbolized  s,  are  supposedly  measured  in  units  such 
that  the  cc  over  a unit  distance  is  0.99.  Suppose,  however,  distances  are  meas- 
ured in  km,  and  symbolized  s'.  Then  the  distance  so  measured,  over  which  the 
cc  is  0.  99,  is  defined  as  the  scale  distance  and  is  symbolized  r (km).  It  is  given 
by 


r = s'/s 


(27) 


Given  two  of  r,  s',  or  s the  third  is  obtained  by  this  equation.  Both  s'  and  r must 
be  measured  in  the  same  units  (for  example,  km),  and  s is  dimensionless. 


1.1  Model  \ with  Seale  Distance  |r) 


Putting  A = 50  in  (22)  the  cc  of  Model  A decreases  as  the  square  of  the  dis- 
tance, thus: 


p(s')  = exp  {-(s'/10r)  } 


(28) 


1.2  Model  B with  Scale  Distance  (r) 

Putting  R = 64  in  (25),  and  using  (24): 


p(  s 


2 f . -1  2 aH  2 

')  - — sin  VI  - it  - a \ 1 - a 

7T  L 


J' 


where 


tr  = s'/128r 


(29) 


1 

I 


1 
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For  a computer  exercise,  the  square  grid  (Figure  2)  was  chosen  to  consist  of 
97  rows  and  97  columns.  To  produce  each  map  random  numbers  (yQ)  were  initially 
assigned  to  each  of  the  97  X 97  grid  points.  The  radius  of  influence  <Rq>  of  Model 
B was  chosen  equal  to  16  units.  The  choices  of  97  X 97  grid  points  and  Rq  - 16 
were  dictated  by  computer  economics. 

Through  (151  and  (20)  for  Model  A or  (23)  for  Model  B,  the  random  number  at 

each  point  (i,  j)  on  the  central  65  rows  and  65  columns  was  replaced  by  y(i,  j). 

The  summations  required  in  (15)  were  extended  only  to  those  797  points  within  the 

16-unit  radius  circle  surrounding  the  point  (i,  j).  For  a radius  of  influence 

R 16  the  cc  in  Model  B between  points  that  are  one  unit  distance  apart  is 
o ’ 

p 0.  96.  To  obtain  the  same  cc  one  unit  distance  apart  in  Model  A the  parameter 
was  made  A 12.  50. 

One  sample  set  of  initial  random  numbers,  processed  first  by  Model  A,  then 
by  Model  B,  produced  the  results  as  shown  (Figure  4).  In  a 9-row  by  9-column 
area,  the  Model  A values  of  y(i,  j)  varied  from  a high  of  2.  4 to  a low  of  -1.  3.  In 
Model  B there  was  a smaller  range  of  y(i,  j)  in  the  same  area,  because  the  spatial 
cc  maintains  a much  higher  value  throughout  the  area.  But  whereas  in  Model  A 
smooth  isopleths  can  usually  be  drawn  around  well-formed  centers  of  low  and  high, 
in  Model  B the  field  is  frequently  erratic.  As  a digression  from  the  main  theme  of 
this  paper,  an  inference  can  be  drawn  from  this  result  with  regard  to  the  spacing 
of  weather  stations.  If  the  process  in  nature  resembles  the  Model  A process  then 


MODEL  A 


p[ 8)«0  074 


MODEL  B 


P(6)*0685 


Figure  4.  Sample  Synoptic  Fields  Produced  By  Model  A 
and  Model  B from  Some  Initial  Set  of  Random  Numbers 
Assigned  to  the  Grid  Points 
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relatively  few  observation  points  will  be  needed  to  complete  the  drawing  of  isopleths 
with  reasonable  accuracy.  If,  however,  the  natural  process  resembles  Model  B, 
then  even  with  closely  spaced  observations  some  significant  events  might  be 
missed  by  the  network  of  stations. 

In  the  example  shown  in  Figure  5 an  initial  set  of  random  numbers  produced  a 
similar  typical  result.  Supposing  that  the  frequency  of  a cloud  cover  corresponds 
to  y < -0.  5,  and  the  frequency  of  rainfall  corresponds  to  y < -1.0,  then  the  Model 
A picture  shows  a large  single  cloud  cover  with  well-defined  boundaries  and  a 
single  large  patch  of  rain.  The  Model  B picture  shows  more  irregularity  in  the 
boundary  of  the  cloud  cover  and  separate  small  patches  of  cloud  and  breaks  in  the 
main  cloud  deck;  also  rainfall  is  represented  by  two  cells  within  the  large  cloud. 

Zf  the  dimension  of  Figure  5 is  of  the  order  of  67  km,  then  Model  B appears  to  be 
a better  simulator  of  the  mesoscale  characteristics  of  clouds  and  rainfall,  espec- 
ially summertime  rainfall. 


MODEL  A 


MODEL  B 


Figure  5.  Comparison  of  Two  Pictures  of  Supposedly 
Cloud  Cover  (single  hatching)  and  Precipitation  (double 
hatching).  The  images  are  stochastically  produced  by 
Model  A and  Model  B from  the  same  initial  set  of  ran- 
dom numbers 


6.  PROBABILITY  OISTRIBI  TION  OF  SPATIAL  MINIMI  M (OR  MAXIM!  M) 

For  a selected  model  with  given  function  F(R)  the  simulation  program  produces 
a horizontal  map  (for  example.  Figures  4 or  5).  Repeated  M times,  it  will  produce 
M different  maps.  At  any  one  point  (i,  j)  of  the  grid  (Figure  2)  the  M values  of 
y(i,  j)  from  the  M maps  will  collectively  have  a Gaussian  distribution,  for  M large. 
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Suppose  (s  + 1)  simultaneous  consecutive  values  of  yU,  3)  along  a single  row 
(j  constant)  are  examined  for  a minimum.  From  M maps  there  will  be  M values 
of  the  (s  * 1)  - station  minimum,  thus  giving  a frequency  distribution  of  the  mini- 
mum along  a line  of  travel  of  length  s units.  Such  simulation  is  used  as  this 
paper's  estimate  of 


P(-Vmin:d)  * 

defined  as  the  probability  that  the  minimum  of  y along  a straight-line  segment  of 

length  s units  will  be  equal  to  or  less  than  ymin. 

Suppose  that  the  (s  + l)2  simultaneous  values  of  a part  of  the  grid,  that  forms 

a square  with  sides  s units  long,  are  examined  for  a minimum.  From  M maps 

there  will  be  M values  of  the  (s  + l)2  - station  minimum,  thus  giving  a frequency 

2 

distribution  of  the  minimum  in  an  area  of  s . Such  simulation  is  used  as  this 
paper's  estimate  of 

P(ymin:s2)  - 

2 

defined  as  the  probability  that  the  minimum  value  of  y throughout  the  area  s will 
be  equal  to  or  less  than  ymin-  The  density  of  stations  is  assumed  great  enough  to 
permit  us  to  consider  the  (s  + l)2  - station  minimum  as  the  minimum  of  continuous 
values  inside  the  area. 

A reversal  of  the  signs  of  all  y in  this  exercise,  yields  the  corresponding 
results  for  the  maximum 

p(ymax:s)  > 

defined  as  the  probability  that  the  maximum  of  y along  a straight-line  segment  of 
length  s will  be  equal  to  or  less  than  ymax.  and 

P(y  ;s2) 

'max 

2 

defined  as  the  probability  that  the  maximum  of  y throughout  the  area  s will  be 
equal  to  or  less  than  ymax- 

It  is  clear  that  a practical  exercise,  as  described  above,  can  produce  only 
approximations  by  simulation  to  the  probabilities  of  spatial  minima  or  maxima. 

The  approximations  will  be  improved  with  a very  fine  mesh  (Figure  2)  or  with  very 
many  rows  and  columns  and  strong  persistence  between  adjacent  points.  However, 
limitations  of  computer  time  made  it  necessary  to  use  relatively  few  rows  and 
columns  for  manageable  sample  sizes,  as  described  below.  For  larger  areas  the 
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sample  size  (M)  was  made  smaller,  and  the  persistence  between  adjacent  points 
was  reduced  by  making  the  parameter  smaller. 


l>.  1 Model  \ Distribution 

With  the  parameter  in  (20)  set  at  A = 50,  then  (22)  gives  the  cc  over  a unit  dis- 
tance (s  = 1)  as  p 0.99.  The  weighting  function  F(R),  beginning  with  1.0  at  (i,  j), 
drops  as  low  as  0.  006  at  a radial  distance  of  16  units  from  point  (i,  j).  For 

17  £ i,  j £ 25  the  maps  were  square,  8 units  to  the  side.  Taking  M = 10,  000,  the 

2 

sample  size  was  large  enough  to  give  values  of  and  ^^min'3  ^ that,  when 

plotted  on  probability  paper,  permitted  the  reasonably  accurate  drawing  of  the  parts 

of  the  curves  of  ymin  in  Figures  6 and  7 for  s £ 8.  The  curves  in  Figures  6 and  7 

are  drawn  for  y . = -3.  5 (0.  5)  3.  0.  3'he  s -scale  on  the  horizontal  axis  is 

min 

logarithmic  and  begins  on  the  left  side  at  s = 1/2.  There  is  also  a uniform  z-scale 
on  the  horizontal  axis,  where 

z = Ins/in2  . (30) 

On  the  ordinate,  ymin  is  considered  to  be  simply  y(0,  1),  that  is,  the  minimum  in 
an  infinitesimal  area  or  on  an  infinitesimally  small  line  segment.  The  cc  is  0.  99 
over  the  distance  s = 1. 

To  extend  the  results  to  larger  areas,  up  to  s - 128  or  z = 7,  the  parameter  in 
(20)  was  set  at  A"  = 12.  5,  and  the  maps  were  made  squares  of  64  units  to  the  side 
(17  £ i,  j £ 81).  With  distances  symbolized  as  s",  the  cc  for  s"  = 1 in  (22)  is 
0.  9608.  This  is  the  same  cc  at  distance  s = 2 when  A = 50.  In  general  s = 2s". 

For  the  part  of  the  curve  where  8 < s £ 128,  the  sample  size  was  made  M = 1000, 
and  the  results  were  plotted  and  drawn  continuous  with  the  previous  segments  for 
s between  1/2  and  8 units. 

To  extend  the  results  to  a still  longer  distance  and  larger  area,  up  to  s = 256 
or  z = 8,  the  parameter  was  made  A"  =3.  125  for  which  s = 4s".  Again  the  sample 
size  was  made  M = 1000,  and  the  results  were  plotted  and  smoothed  with  the  pre- 
vious results. 

In  Figures  6 and  7 the  y(0,  l)-scale  is  the  scale  of  ymjn>  and  the  values  for 

P(y  . ;s)  or  P(y  . ;s 2>  are  read  on  the  left-hand  scale.  If  the  signs  of  y(0,  1)  are 
m in  m in  » 

reversed  then  the  values  of  P(y  ;s)  or  P(y  ;s  j are  read  on  the  right-hand 

max  max 

vertical  scale. 

6.2  Model  R Spatial  Minimum  (or  Maximum) 

In  this  model  it  would  not  be  feasible  to  have  a radius  of  influence  R = 64, 

o ’ 

since  it  would  require  close  to  13,  000  random  numbers  initially  to  obtain  a value  of 
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The  other  is  the  uniform  scale  of  z 
;s  are  for  values  of  y • - -3.5  (0. 


v at  each  point  (i,  j).  With  the  radius  of  influence  set  at  R " = 16,  then  (24)  gives 
the  cc  over  a unit  distance  (s"  1)  equal  to  0.  9602.  This  is  the  same  cc  as  at 

distance  s 4 when  R(  64.  For  s 1,  the  cc  remains  0.99.  In  general  s 4s". 
The  weighting  function  is  1.  0 throughout  the  circular  area  of  radius  Ro".  For 
17  -s  i,  j £ 25  the  maps  were  square,  8 units  to  the  side.  Taking  M 75,  000,  the 
values  of  and  P(>min’s^^  were  estimated  to  a high  degree  of  accuracy, 

which  permitted  drawing  the  curves  for  s <32  (Figures  8 and  9). 

To  extend  the  results  to  larger  areas,  up  to  s 256  or  z 8,  the  sample  size 
had  to  be  reduced  (M  = 5000)  to  obtain  estimates  within  a reasonable  time  of  com- 
puter operation.  The  synoptic  maps  were  squares,  64  units  to  the  side  for  17  £ i, 

) £ 81.  The  results  for  the  minima  were  plotted  and  smoothed  with  the  previous 
results . 

The  scales  in  Figures  8 and  9 have  similar  descriptions  as  those  in  Figures 
6 and  7 and  the  charts  are  used  in  similar  fashion. 

ti.-i  Model  It  lliehe-i  Minimum  lor  Lowest  Maximum!  in  a Fractional  Loverage 
of  an  \rea 

Each  model,  as  described  so  far,  provides  an  estimate  of  the  probability  of  a 
minimum  or  maximum  condition  in  a region  of  stated  size.  But  there  is  still  a 
requirement  to  estimate  the  probability  that  a certain  percentage  or  fraction  of  a 
given  area  will  experience  a critical  event  or  a threshold  condition.  For  example, 
in  the  case  of  heavy  24-hr  rainfall  (>S0  mm),  it  might  be  desirable  to  know  the 

probability  that  one-fourth  or  one-half  of  a city  will  receive  this  amount,  compared 

2 

with  the  lesser  probability  that  the  whole  city  of,  say,  100  km  will  be  drenched. 

2 

In  a horizontal  map  (for  example,  Figure  4)  suppose  that  the  (s  + 1)  simul- 
taneous values  (y)  are  sorted  from  the  largest  to  the  smallest  (y'h  Then  the  first 

f(s  + 1)“  values  of  y'  will  be  equal  to  or  greater  than  their  last  value.  Any  other 

2 

selection  of  a fraction  (f)  from  the  (s  + 1)  values  of  y will  have  a lower  minimum. 

2 

From  M maps  there  will  be  M values  of  y'{f(s  + 1)  },  to  provide  a frequency  dis- 

2 

tribution  of  the  maximized  minimum  in  the  fractional  coverage  (f)  of  an  area  (s  ). 
Such  simulation  is  used  as  this  paper's  estimate  of 

P(ymaxmin:s2' f)  - 

2 

defined  as  the  probability  that,  in  a fraction  f of  the  area  s , the  highest  possible 
minimum  value  of  y is  y„  . . A reversal  of  the  signs  of  all  y's  in  this  exercise 
yields  the  corresponding  results  for  the  minimized  maximum 


P(y  • ;s  , f)  , 
•'minmax  ’ ’ 


Model  B Probability  Estimates  of  the  Mini-  Figure  9.  Model  B Probability  Estimates  of  the  Mini 
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defined  as  the  probability  that,  in  a fraction  f of  the  area  s2,  the  lowest  possible 
maximum  value  of  y is  y . _ 

Model  B was  chosen  for  further  exploration  by  Monte  Carlo  procedure,  since 
it  proved  to  be  the  more  effective  model  on  precipitation  (see  Section  7).  In  each 
synoptic  map,  consider  the  smallest  square  in  the  grid  (any  square).  With  four 
values  of  y,  one  at  each  corner  of  the  smallest  square,  when  ordered  by  size,  the 
ith  value  can  be  considered  to  be  the  ith  largest  (i  1,  4)  in  a sample  of  four  values 
taken  from  the  infinite  population  in  the  unit  area.  Consider  the  square  of  2 units 
to  each  side.  When  the  9 associated  values  are  ordered  by  size,  the  ith  value  is 

considered  to  be  the  ith  highest  in  a sample  of  9 taken  from  the  infinite  population 

2 

of  the  area.  In  general,  in  a square  of  s units  to  one  side,  when  the  (s  + 1)  values 

2 

are  ordered  by  size,  the  ith  value  is  the  ith  highest  in  a sample  of  (s  +•  1)  values. 

7 

The  fraction  (f)  of  the  area  represented  by  the  i values  is  estimated  by 

F = (i  - 3/8)/  { (s  + l)2  + 1/4}  (31) 

In  the  Monte  Carlo  exercise  the  number  of  synoptic  maps  generated  from  initially 

random  numbers  was  M = 5000.  Each  man  was  64  units  to  one  side,  so  that  seven 
2 2 2 2 2 ° 2 

square  areas,  of  sizes  1,2,4,8,13,  32“,  64  could  be  studied.  But  the 
field  was  produced  for  a radius  of  influence  R ' = 16,  so  that  the  cc  between  sta- 
tions one  unit  apart  (s'  = 1)  was  0.9602.  That  is,  the  smallest  square  is  actually 
4 units  to  a side  (s  = 4),  when  measured  by  that  element  of  length  over  which  the 
cc  is  0.  99. 

Consider  an  elementary  area  of  one  square  (s'  = 1 or  s = 4)  with  the  four  cor- 
ner values  of  y ordered  by  size.  For  M = 5000  the  cumulative  distribution  of 

values  of  y for  the  first,  second,  third,  and  fourth  values  were  plotted 

(Figure  10)  at  fraction  (f)  of  the  area  given  by  (31).  On  the  left-hand  ordinate 
(Figure  10)  the  probability  is  that  of  the  highest  minimum  in  an  infinitesimal  frac- 
tion of  the  area,  in  other  words  the  areal  maximum,  as  obtained  from  Figure  9 
(at  s = 4 or  z = 2).  On  the  right-hand  side  the  cumulative  probability  is  that  of  the 
minimum  of  the  whole  area,  again  obtained  from  Figure  9 (at  s = 4 or  z = 2).  Thus, 
each  y-curve,  y = -3.  5 (0.  5)  3.  5,  was  drawn  through  six  points,  at  each  end  and  at 
f = 0.  147,  0.  382,  0.  618  and  0.  853. 

Consider  the  larger  area  (s'  = 2 or  s = 8)  of  4 squares  composed  of  9 grid 
values  of  y that  are  then  ordered  by  size.  In  the  large  sample,  M = 5000,  the 
distributions  of  each  of  the  9 values  of  y'  were  plotted,  in  accordance  with  (31)  at 
f = 0.  0757,  0.  197,  • • •,  and  0.  932,  along  with  the  distribution  of  the  minimum 


Blom,  Gunnar  (1958)  Statistical  Estimates  and  Transformed  Beta-Variables, 
John  Wiley  and  Son=  New  York. 
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on  the  right-hand  aide  and  the  maximum  on  the  left-hand  aide,  aa  found  in  f igure  9 

at  a 8 (or  z 3).  The  chart  ao  produced  waa  similar  to  Figure  10.  Likewise, 

2 2 2 2 2 

similar  charts  were  drawn  for  the  areas  4,8,  16  , 32“,  and  64".  The  isopleths 
of  all  seven  such  charts  provided  points  for  the  plotting  of  Figures  11  to  19  for  the 
highest  1/10,  2/10,  •••,  highest  9/10.  The  last  chart  of  this  series,  for  the 
highest  10/10,  would  be  identical  with  Figure  9. 

Figures  11  to  19,  together  with  Figure  9,  comprise  the  substance  of  this 
section  and  provide  the  means  for  determining  the  probabilities  of  fractional  cover 
of  a given  area  by  a describable  event.  At  first  glance  the  charts  seem  to  present 
a few  surprises,  but  upon  further  reflection  the  results  are  reasonable.  For  exam- 
ple, an  event  that  has  a 50  percent  chance  of  occurring  at  one  single  point  will 
occur  with  greater  probability  (read  right-hand  scale)  over  1/10  to  4/10  of  the 
area.  On  the  other  hand,  this  probability  must  decrease  for  fractions  of  the  area 
greater  than  une-half.  Again,  for  a rare  event  of,  say,  1 percent  single-point 
probability,  the  chances  are  improved,  to  as  much  as  3 percent  of  its  occurring  in 
1/10  of  the  area  — up  to  a certain  optimum  areal  size.  But  as  the  area  is  further 
enlarged,  the  probability  of  the  rare  event  covering  as  much  as  1/10  of  the  increas- 
ingly larger  area  begins  to  vanish. 

ippLirATiov* 

For  the  most  direct  application  of  either  Model  A or  Model  B to  a weather 
element  (X)  the  input  consists  of  the  single-point  frequency  distribution,  P(x), 
plus  one  additional  parameter:  the  scale  distance  (r).  Through  the  graphs 
(Figures  6 through  9,  11  through  19)  the  information  is  used  to  estimate  the  prob- 
ability of  occurrence  along  a line  of  given  length,  or  within  an  area  of  given  size 
or  a fraction  of  the  area. 

For  a one-to-one  correspondence  between  x and  y,  as  defined  in  (1),  let  us 
further  define: 

P = P(x  . ) = P(y  . ),  the  probability  that  x will  not  exceed  x . at  a single 
o min  •'min  ’ 1 J min  6 

selected  point. 

P P(x  . ;s')  = P(y  . ;s),  the  probability  that  the  minimum  of  x along  a line 
s min  •'min  ’ r J 6 

of  length  s'  (km)  will  not  exceed  x . . 

2 2 min 

P P(x  . ;s'  ) - P(y  . ;sM,  the  probability  that  the  minimum  of  x through- 
a mm  •'min  2 2 

out  the  area  s1^  (kin  ) will  not  exceed  x 
2 2 ^ ^ 

Pf  = P(xmaxmin:s'  * f)  = P<ymaxmin;s  ’ f)>  the  Pr°bability  that  the  maximized 

minimum  of  x in  a fraction  f of  the 
2 2 

area  s'  (km  ) will  not  exceed 
xmaxmin’ 
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Probability  Estimates  of  the  High-  Figure  15.  Model  B Probability  Estimates  gf  the  High 


( 01/8  ) 


Probability  Estimates  of  the  High-  Figure  19.  Model  B Probability 


When  the  subject  is  spatial  maximum  we  define: 

P P(x  ) = P(y  ),  the  probability  that  x will  not  exceed  x at  a sin- 

o max  •'max  ’ J max 

gle  point,  essentially  the  same  definition  when 

P = P(x  . ). 
o min 

P = P(x  ;s')  - P(y  ;s),  the  probability  that  the  maximum  of  x along  a 

s max  max 

line  of  length  s'  (km)  will  not  exceed  x . 

2 9 max 

P = p(x  ; s '*)  = P(y  ;s  ),  the  probability  that  the  maximum  of  x through- 

a max  max  2 2 

out  the  area  s'  (km  ) will  not  exceed  x . 

2 2 ax 

P.  = P(x  . ;s'  ,f)  = P(y  . ;s  , f),  the  probability  that  the  minimized 

f minmax  ’ 'minmax  ’ ’ 1 J 

maximum  of  x in  a fraction  f of  the 
2 2 

area  s'  (km  ) will  not  exceed 


In  the  following  examples,  procedures  were  followed  to  produce  an  estimate 
of  scale  distance  (r),  which  then  becomes  a telling  characteristic  of  the  weather 
element  under  investigation.  In  Model  A or  Model  B the  larger  the  value  of  r the 
greater  is  the  horizontal  persistence  of  the  element.  The  first  example,  below,  is 
of  the  macroscale  100-mb  wintertime  temperature  field  in  the  northern  hemisphere. 
The  second  example  is  of  the  horizontal  extent  of  New  England  rainfall.  With 
respect  to  rainfall  a difficulty  arises  because  of  its  integration  over  a time  inter- 
val, like  one  hour  or  one  full  day.  The  third  example,  on  the  other  hand,  consist- 
ing of  a study  of  PPI  radar-scope  pictures,  is  considered  to  be  a more  logical 
application  of  the  model  of  this  paper,  since  each  picture  is  closer  to  an  instan- 
taneous representation  of  the  horizontal  spread  of  the  meteorological  condition. 

Suppose,  for  an  element  (X)  that  the  climatic  frequency  is  used  as  an  estimate 
of  the  probability  P . Suppose  additionally  that  the  frequency  with  which  x is  the 
minimum  along  a line  of  length  s'  is  used  to  estimate  probability  P . Each  value 

of  P defines  a value  of  y . on  the  ordinate  scale  of  Figure  6 when  Model  A is 
o min 

used,  and  Figure  8 when  Model  B is  used.  The  value  of  s is  found  where  the  ymin 
curve  (generally  interpolated)  intersects  the  probability  P . More  practically  the 
value  of  z (Eq.  30)  is  found  on  the  uniform  scale  of  the  horizontal  axis,  and  then  s 
is  obtained  through  the  transformation 


Likewise  if  the  frequency  is  P&,  with  which  is  the  minimum  throughout  the 

area  s'2,  the  value  of  s is  found  in  Figure  7 when  Model  A is  used  and  in  Figure  9 

when  Model  B is  used,  where  the  y . curve  intersects  the  probability  P . Again, 

m in  a. 

if  the  frequency  is  Pf,  with  which  xmaxmjn  is  the  highest  minimum  in  a fraction  f of 
the  area  s'2  for  f = 1/10,  • • • , 9/10,  the  value  of  s or  z in  Model  B is  found  in 
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Figure  11.  Figure  12,  •••,  or  Figure  19,  respectively,  where  the  ymaxrnin  curve 
intersects  the  probability  P^..  Each  ratio  of  s'  to  s is  an  estimate  of  the  scale  dis- 
tance  (r). 

For  m(a)  values  of  X in  the  ath  areal  size,  and  for  n areal  sizes,  there  will  be 
N estimates  of  the  scale  distance  (r)  where 


n 

N = m(a> 

a = l 


(33) 


The  proximity,  to  each  other,  of  all  individual  estimates  of  r is  a measure  of  the 
effectiveness  of  the  model.  The  order  of  magnitude  of  the  values,  as  in  the  exam- 
ples below,  could  appear  to  reveal  a damaging  lack  of  similarity  in  the  judgment 
of  the  reader.  But  the  models  were  useful  in  application  and  a single  "best"  esti- 
mate was  well  justified. 

How  should  the  "best"  single  estimate  of  r be  obtained?  Since  r is  the  ratio 
of  s'  to  s,  the  errors  on  the  low  side  are  bounded  by  zero,  but  on  the  high  side  the 
errors  are  unbounded.  On  the  other  hand,  In  r is  equally  unbounded  on  both  the 
low  and  the  high  side.  Moreover,  since  the  graphs  (Figures  5 through  9,  11  through 
19)  are  constructed  with  a uniform  scale  of  z,  which  is  proportional  to  in  s,  the 
errors  in  estimating  Ins  tend  to  be  more  equally  divided  above  and  below  a true 
value  than  they  would  be  in  estimating  s.  Consequently,  it  is  preferable  to  take  the 
average  oflnrfdns1  - in  s)  than  the  average  of  r per  se.  This  is  the  same  as 
taking  the  geometric  mean  of  r,  as  follows: 


m(a),  n 

n (»a'/»al) 

i=l,  a=l 


1/N 


(34) 


In  practice  each  value  of  z is  obtained  from  the  appropriate  curve  (Figures  6 
through  9,  11  through  19)  by  using  each  pair  of  PQ  and  Pg,  or  Pft,  or  Pf.  For  the 
ath  area  the  average  value  of  z is  given  by 


m(a) 

z - £ Zj/m(a)  , (35) 

1=1 

wnich  yields 

s(a)  = 2Z  , (36) 
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which,  when  divided  into  s'(a),  yields 


1 


r(a)  ; s'(a)/s(a) 


(37) 


Then,  lastly,  a single  estimate  of  r is  obtained  by 


n 

{r(a)}m(a) 

a 1 


(38) 


7.1  Example  1 100-mb  Route  Temperatures 

Climatological  maps  of  the  northern  hemisphere  temperature  at  supersonic 

g 

transport  altitudes  have  been  published  by  season  in  percentiles  from  2 to  98. 

Each  map  is  accompanied  by  a graph  that  provides  an  estimate  of  the  probability 
that  the  temperatures  on  a given  route  and  midlatitude  will  remain  above  the  tem- 
peratures that  are  on  the  map  along  that  route.  For  the  100-mb  level,  winter- 
time, these  probabilities  (P^)  are  as  shown  in  Table  2 for  routes  of  varying  length 
from  1000  km  to  9000  km,  with  midlatitude  45°N.  The  probability  estimates  (Pg) 
had  been  made  by  actual  sampling  of  the  data  on  synoptic  charts. 

Model  A (Figure  6),  when  applied  to  the  values  of  PQ  and  Pg  (Table  2),  yields 
values  for  z,  one  value  corresponding  to  each  pair  of  PQ  and  Pg.  Consequently, 
through  (32),  (27),  (37),  Model  A yields  an  estimate  of  scale  distance,  r - 289  km, 
which,  when  applied  in  (28),  gives  a decrease  of  cc  with  distance  as  shown 

9 

(Figure  20).  Also  shown  is  Bertoni  and  Lund's  curve  for  the  decrease  of  temper- 
ature cc  between  pairs  of  stations  at  the  16-km  level  (near  100-mb).  In  contrast 
with  Bertoni  and  Lund's  data  this  author's  data  were  taken  from  smoothed  synoptic 
analyses  of  100-mb  surfaces,  which  would  tend  to  increase  the  apparent  horizontal 
relationship  between  stations.  Moreover,  the  Bertoni  and  Lund  temperatures  were 
derived  from  radiosonde  data.  Hence  the  temperatures  on  constant -level  surfaces 
(rather  than  constant-pressure)  might  have  additional  errors,  causing  faster  cc 
decay  and  calling  for  a smaller  scale  distance  (see  below). 

Model  B (Figure  8)  yielded  the  other  estimates  for  z as  shown  in  Table  2 and 
a consequent  scale  distance  r = 85.  6 km.  When  used  in  (29)  and  (24)  the  decay  of 
cc  becomes  as  shown  (Figure  20).  Clearly  Model  A is  preferable  for  the  macro- 
scale upper-air  temperatures. 

8.  Gringorten,  1. 1.  and  Tattelman,  P.  (1970)  Point  and  Route  Temperatures  for 

Supersonic  Aircraft,  AFSG,  No.  233,  AFCRL-70-0420,  Bedford,  MA. 

9.  Bertoni,  E.A.  and  Lund,  I.A.  (1964)  Winter  Space  Correlations  of  Pressure. 

Temperature  and  Density  to  16  km,  ERP,  No.  75,  AFCRL-64-1020, 

Bedford,  MA. 


36 


I 


Table  2.  The  Cumulative  Probability  (P^)  of  the  Minimum  Temperatures 
(percentile  PQ)  at  100  mb.  Wintertime,  Along  a Route  Whose  Midlatitude 
is  45°N.  The  probability  that  the  temperatures  enroute  will  equal  or 
exceed  the  percentile  ( P0)  is  (1  - Ps).  For  the  last  three  columns  see 
text 


Percentile 

P 

o 

. 

Temps  on 
map  (°C) 
Equator 
to 

N-pole 

Route 
length 
(km) 
s * 

P 

S 

Model  A 

z 

Mode  1 B 

When 
r=185  km 
P 

s 

0.02 

-88 

1481 

0.03 

1.62 

3.16 

0.045 

to 

2130 

0.04 

2.61 

4.26 

0.055 

-55 

3148 

0.05 

3.21 

4.90 

0.069 

4908 

0.07 

4.17 

5.78 

0.089 

7223 

0.  10 

4.97 

6.65 

0.109 

8519 

0.11 

5.18 

6.81 

0.128 

0. 10 

-82 

1074 

0.12 

1.82 

2.25 

0.170 

to 

1824 

0.15 

2.05 

3.71 

0.205 

-53 

3148 

0.20 

3.22 

5.00 

0.268 

4815 

0.25 

3.91 

5.83 

0.336 

6112 

0.30 

4.44 

6.33 

0.380 

7130 

0.35 

4.80 

6.74 

0.412 

8704 

0.40 

5.16 

7.06 

0.470 

0.25 

-81 

1 1 30 

0.  30 

1 . 10 

2.5 

0.375 

to 

1611 

0.35 

2.22 

3.70 

0.420 

-50 

2797 

0.40 

2.93 

4.58 

0.50 

3704 

0.50 

3.90 

5.73 

0.565 

4630 

0.55 

4.21 

6.09 

0.620 

6112 

0.60 

4.53 

6.44 

0.695 

7408 

0.70 

5.08 

7.10 

0.745 

9075 

0.80 

5.63 

7.65 

0.800. 

O 

cn 

O 

-80 

1296 

0.60 

2.00 

3.51 

0.66 

to 

1778 

0.65 

2.68 

4.37 

0.707 

-48 

2426 

0.70 

3.17 

5.05 

0.76 

3704 

0.80 

4.01 

6.00 

0.84 

6297 

0.90 

4.82 

6.85 

0.927 

8704 

0.95 

5.36 

7.36 

0.964 

i r 
C" 

o 

-79 

982 

0.80 

1.26 

2.80 

0.85 

to 

1463 

0.85 

2.40 

4.30 

0.884 

-45 

0.90 

3.24 

5.26 

0.934 

4630 

0.95 

4.08 

6.18 

0.975 

7778 

0.99 

5.14 

7.15 

0.9945 

o 

O' 

o’ 

-78 

1306 

0.95 

2.28 

3.88 

0.965 

to 

3148 

0.98 

3.47 

5.38 

0.9917 

-43 

4352 

0.99 

4.00 

5.93 

0.9967 

0.98 

-77  to 

1389 

0.99 

1.76 

3.35 

0.9953 

-41 

■ 
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Figure  20.  Decrease  of  16-km  Temperature  cc 
with  Distance  Between  Stations,  as  Found  by 
Bertoni  and  Lund  (1964),  and  the  Model  A and 
Model  B Estimates  of  the  Decrease  on  the  100- 
mb  Surface 


An  application  of  Figure  6 or  Figure  7 together  with  an  assigned  value  of  the 
scale  distance  (r)  provides  a method  for  obtaining  the  probability  of  the  minimum 
percentile  of  route  temperature,  or  the  minimum  in  an  area  of  given  size.  It  is 
still  necessary  to  have  geographic  charts  showing  w'hat  temperatures  correspond 
to  the  percentiles  of  0.  02  • • • , 0.  99. 

Suppose  the  problem  is  to  find  versus  s',  knowing  the  cc  decrease  with 
distance  as  given  by  Bertoni  and  Lund  (Figure  20).  For  Model  A the  Bertoni  and 
Lund  curve  yields  a geometric  mean  value  of  the  scale  distance,  r = 185  km, 
which  is  36  percent  shorter  than  the  estimate  using  the  100-mb  data.  Hence,  with 
less  horizontal  persistence,  it  is  to  be  expected  that  low  minima  will  be  predicted 
with  greater  probability.  Using  r 185  km  to  obtain  values  of  z from  (27)  and 
(30)  to  correspond  to  the  values  of  s'  in  Table  2,  and  entering  z in  Figure  6,  esti- 
mates (P  .)  were  obtained  as  shown  in  the  last  column  of  Table  2.  On  the  average, 
these  latter  probability  estimates  are  13  percent  higher. 

7.2  Example  2:  Rainfall  frequencies 

Recent  monthly  publications  of  the  National  Weather  Service  reveal  that  there 
are  more  than  3 50  stations  at  which  24 -hr  rainfall  is  observed,  once  per  day, 
throughout  New  England.  Such  a network  does  provide  a good  areal  coverage  of 
precipitation.  But  it  confronts  this  study  with  several  difficulties.  First  of  all  the 
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cooperative  observers  do  not  have  a single  time  at  which  to  read  the  24 -hr  cumula- 
tion of  precipitation.  At  most  stations  the  rain  gauge  is  read  at  0800  Local  Time. 
But  many  observers  read  the  amount  in  the  afternoon,  some  at  midnight.  Hence, 
a survey  of  rainfall  amounts,  to  obtain  the  24 -hr  maximum,  would  be  faulty  unless 
the  data  were  used  only  from  stations  taking  observations  at  the  same  time. 

Some  100  stations  report  precipitation  amounts  hourly  and  they  provide  24-hr 
amounts  ending  at  midnight.  With  the  study  limited  to  these  stations,  the  maximum 
daily  precipitation  was  recorded  in  each  of  the  areas  shown  in  Table  3. 


Table  3.  Showing  Areas  of  Increasing  Size,  Beginning  with  Massachusetts 
and  Ending  with  All  of  New  England.  The  distance  s' (a)  is  the  square  root 
of  the  ath  area 


(1)  New  England  January  24-hr  Precipitation 

The  precipitation  data  presented  in  the  National  Environmental  Data  Service 
monthly  bulletins  (U.  S.  Dept,  of  Commerce,  1952-1972)  were  surveyed  to  obtain 
the  climatic  frequency  (P  ) of  24-hr  precipitation  and  the  frequency  of  maximum 

(P  , for  a = 1,  2,  3,  4)  in  the  successively  larger  areas  of  New  England  beginning 
a 2 

with  the  21,  358  km  of  Massachusetts  (Table  4).  The  true  climatic  frequencies  of 
precipitation  are  assumed  to  be  the  same  for  all  four  areas  and  therefore  the  dif- 
ferences in  Pq  shown  in  Table  4 are  assumed  to  be  sampling  errors  from  true 
probabilities.  For  the  intended  purpose  of  finding  the  parameter  (r),  each  fre- 
quency (P  ) is  used  in  conjunction  with  the  frequency  of  the  maximum  precipitation 
(P  ) for  its  area.  Table  4 presents  only  a fraction  of  the  cumulative  frequencies 
used  for  24-hr  rainfall.  Actually  used  were  the  cumulative  frequencies  as  follows: 
none  or  trace,  0.  25(0.25)2.3,  3.  6,  4.8,  7.4(2.  54)27.  7,  37.  8(12.  70)88.  6 mm. 

Each  value  of  PQ  was  entered  on  the  left-hand  ordinate  of  Figure  9,  although 
read  on  the  right-hand  scale.  The  Pq  curve,  usually  interpolated,  was  followed 
until  the  value  of  P was  reached,  then  the  value  of  z on  the  abscissa  was  read 
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States 

Area 

no. 

(a) 

Area 

km2 

s ' (a) 
km 

MA 

i 

21358 

146.1 

MA,RI,CT 

2 

37455 

193.5 

MA,RI,CT,NH,VT 

3 

86377 

293.9 

MA,RI,CT,NH,VT,ME 

4 

172294 

415.1 

Table  4.  A Partial  Summary  of  the  Cumulative  Relative  Frequency  and  the  Resulting  Value  of  z (see  text) 
of  New  England  January  24-hr  Precipitation  in  21  Years  (1952-1972).  The  frequencies  (P0)  were  based  on 
the  records  of  those  stations  and  years  when  the  31-day  records  were  complete.  The  symbols  used  are  as 
defined  in  the  text.  Areas  1 to  4 are  defined  in  Table  3 


and  recorded  (Table  4).  For  the  greater  extremes  of  precipitation  the  observed 
frequencies  became  erratic.  Hence  the  values  of  were  obtained  only  up  to  the 
frequency  of  50.  a mm  making  m(a)  23  for  a 1,  2,  3,  or  4.  Through  (35),  (3  0), 
(37)  the  final  value  of  the  parameter  in  Model  B for  this  region,  month,  and  24 -hr 
precipitation  was  estimated  at  r 9.  45  km. 

As  stated  above,  the  value  of  r 9.45  km  is  the  estimate  of  the  scale  distance, 
over  which  the  cc  is  0.  99.  Since  Model  B was  used,  the  cc  decreases  with  distance 
as  given  by  (24)  and  (29),  as  shown  in  Figure  21.  When  Model  A (Figure  7)  was 
used  in  a parallel  exercise,  the  scale  distance  r was  30.35  km,  which  when  used 
with  (28)  produced  the  cc  decay  shown  by  the  second  curve  in  Figure  21.  With 
Model  A there  is  a rapid  decrease  of  the  cc  with  distance  between  stations,  after  a 
separation  of  100  km.  To  demonstrate  that  this  is  an  unacceptable  weakness  in 
applying  Model  A,  rather  than  Model  B,  the  precipitation  data  were  surveyed  by 
pairing  the  stations  shown  in  Table  5 to  obtain  cc  more  directly. 


Figure  21.  Decrease  of  cc  with  Distance  (s'  km) 
Between  Precipitation  at  Stations  in  New  England 
in  January.  The  X's  are  plotted  at  six  sample 
estimates  of  tetrachoric  cc  between  the  precipi- 
tation of  five  pairs  of  stations,  based  on  1952- 
1972  data 


For  each  pair  of  stations  there  was  a joint  frequency  of  24-hr  precipitation 
equal  to  or  less  than:  trace,  2.3  mm,  4.8  mm,  7.4  mm,  9.9  mm  and  12.4  mm, 
as  well  as  the  unconditional  frequency  of  these  amounts  at  either  station.  An 
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Table  5.  The  Pairs  of  Stations  that  Were  Selected  to  Find  the  cc  Between 
Their  24-hr  Precipitation  Ending  at  Midnight.  The  January  cc's  are 
shown  plotted  in  Figure  21,  the  July  cc's  are  plotted  in  Figure  22 


First  Station 

Second  Station 

Distance  apart 

Worcester,  MA 

Amherst,  MA 

56  km 

Boston,  MA 

Hartford,  CT 

144  km 

Blue  Hill,  MA 

Augusta,  ME 

256  km 

Hartford,  CT 

Augusta,  ME 

355  km 

Hartford,  CT 

Caribou,  ME 

678  km 

application  of  the  bivariate  normal  distribution10,  11  yielded  six  corresponding 
estimates  of  the  cc  for  each  pair  of  stations,  as  shown  plotted  in  Figure  21.  (The 
cc  obtained  in  this  manner  is  sometimes  known  as  tetrachoric  cc.) 

Clearly  Model  A,  when  the  scale  distance  is  determined  as  above,  does  not 
anticipate  the  high  cc's  between  stations  more  than  100  km  apart.  Model  B also 
tends  to  give  underestimates,  although  they  are  acceptable  for  estimating  order  of 
magnitude  of  the  cc  as  a function  of  distance  of  separation. 

(2)  New  England  July  24-hr  Rainfall 

The  same  source  of  data  was  used  for  the  survey  of  the  frequency  distribution 
of  July  rainfall  and  its  areal  maximum.  Table  G is  a partial  list  of  the  results. 
The  frequency  distribution  of  single-station  rainfall  (P  ),  in  this  case,  is  an  aver- 
age for  some  12  widely  scattered  stations  in  New  England.  The  values  of  z in 
Table  G were  obtained,  through  use  of  Figure  9,  up  to  the  frequency  of  63.  2 mm, 
making  m(a)  = 9 for  a = 1,  2,  3,  or  4.  The  value  of  s'(a)  is  the  square  root  of  the 
ath  area.  Through  application  of  (35),  (36),  (37),  and  (38)  the  final  value  of  the 
scale  distance  for  New  England  July  24-hr  rainfall  was  estimated  at  r 5.  19  km 
with  Model  B. 


10.  National  Bureau  of  Standards  (1959)  Table  of  the  Bivariate  Normal  Distribution 

Function  and  Related  Functions.  NBS  Appl.  Math.  Series,  No.  50,  Govt. 
Printing  Office,  Washington,  D.  C. 

11.  Gringorten,  1. 1.  (1971)  Modelling  conditional  probability,  J . Appl.  Meteorol. , 

10:646-657. 
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When  r 5.  19  km  with  Model  B the  cc  decreases  with  distance  as  shown 
(Figure  22).  When  Model  A was  used  in  a parallel  exercise  the  scale  distance 
was  estimated  at  r 17.  63  km  and  resulted  in  the  second  curve  for  cc  decrease 
with  distance  (Figure  22). 


DISTANCE  ( S'  km) 


Figure  22.  Decrease  of  cc  with  Distance  (s'  km) 

Between  Rainfall  at  Stations  in  New  England  in 
July.  The  X's  are  plotted  at  five  sample  estimates 
of  tetrachoric  cc  between  the  rainfall  at  four  pairs 
of  stations,  based  on  1952-1972  data 

Shown  also  in  Figure  22  are  the  cc's,  five  for  each  separation  of  stations,  as 
determined  by  pairing  stations  to  obtain  the  joint  frequencies  of  July  24-hr  rainfall, 
dichotomously  dividing  their'distributions  at  a trace,  2.3  mm,  4.8  mm,  9.9  mm, 
and  25.  1 mm  and  applying  the  bivariate  normal  distribution  to  obtain  tetrachoric 
cc. 

Again  Model  B tends  to  underestimate  the  cc‘s  slightly,  while  Model  A is 
unacceptable.  For  the  two  stations,  Hartford  and  Caribou,  separated  by  678  km, 
four  of  the  cc's  were  estimated  at  zero  (or  less). 

(3)  New  Englard  January  1-hr  Precipitation 

Since  the  above  exercises  were  conducted  with  24-hr  precipitation,  it  is  pos- 
sible that  the  time -averaging  effect  is  to  increase  the  cc  between  stations,  which 
would  also  imply  an  increase  in  the  scale  distance  (r).  A rainfall  pattern  at  one 
station  frequently  becomes  the  pattern  later  in  the  day  at  a station  east  of  the  first 
station. 
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To  examine  the  impact  of  duration  on  the  results,  a more  laborious  project 
was  undertaken  with  the  hourly  data.  Offhand  the  data  base  and  data  handling 
would  seem  to  be  24  times  as  great,  but  since  the  number  of  no-rain  cases  is  pro- 
portionately greater  among  the  1-hr  statistics,  the  work  was  still  manageable  with 
undergraduate  help  from  Kegis  College,  Weston,  MA. 

The  Hourly  Precipitation  Data  were  surveyed  for  the  frequency  of  the  maximum 
one-hour  precipitation  for  all  24  hours  of  the  day,  disregarding  any  diurnal  cycle 
that  might  exist  in  the  New  England  rainfall.  The  results  for  the  four  successively 
larger  areas  (Table  3)  are  shown  in  Table  7.  The  unconditional  cumulative  fre- 
quencies (P  ) were  obtained  from  the  frequencies  at  five  widely  scattered  stations 
(Boston,  MA ; Burlington,  Vermont;  Caribou,  Maine;  Hartford,  Conn.  ; Pittsfield, 
MA)  that  had  complete  records  and  were  considered  "typical"  of  the  New  England 
rainfall  patterns.  The  values  of  z were  obtained  up  to  the  cumulative  frequency  of 
9.91  mm,  making  m(a)  14.  The  resulting  estimate  of  the  scale  distance,  for 
New  England  January  1-hr  precipitation  was  r - 5.02  km.  This  value  inserted  in 
(29)  and  (24)  would  give  a cc  betwmen  the  precipitation  at  Hartford,  Conn,  and 
Portland,  Maine,  distance  s’  = 286  km  apart,  the  value  of  p(s')  0.45.  By  the 

tetrachoric  method  for  dichotomously  classified  hourly  rainfall  the  cc  varied  be- 
tween 0.  31  and  0.  72  with  an  arithmetic  average  of  0.  58.  Again  Model  B tended  to 
underestimate  the  cc,  but  provided  a proper  order  of  magnitude. 

(4)  New  England  July  1-hr  Rainfall 

The  exercises  on  24-hr  rainfall  supported  a preconception  that  the  rainfall 
patterns  are  more  widespread  and  hence  more  persistent  horizontally  in  midwinter 
than  in  midsummer.  The  24-hr  scale  distance  was  9.45  km  in  January  compared 
with  5.  19  km  in  July.  To  see  if  the  same  is  true  with  the  shorter  duration  in  rain- 
fall, the  exercise  on  hourly  rainfall  in  New  England  was  repeated  using  July  records 
for  the  same  21  years  (1952-1972),  again  disregarding  any  diurnal  cycle  (which 
probably  does  exist  in  summer).  The  results  are  shown  in  Table  8.  But  because 
intensive  rainfall  is  much  more  frequent  in  summer  than  in  winter,  the  values  of  z 
were  obtained  up  to  the  cumulative  frequency  of  27.  7 mm,  making  m(a)  = 21.  The 
resulting  estimate  of  the  scale  distance  for  New  England  July  1-hr  rainfall  was 
r = 3.  61  km  which  again  supports  the  view  that  the  horizontal  pattern  is  less  per- 
sistent in  summer  than  in  winter.  But,  winter  or  summer,  1-hr  rainfall  patterns 
are  expectedly  less  widespread  than  the  24-hr  patterns,  as  revealed  by  the  param- 
eter sizes. 

7.3  Example  3:  PPI-Scope  Radar  Statistics 

It  is  desirable  to  work  on  synoptic  situations  that  essentially  eliminate  time  as 
a factor,  so  that  the  study  can  be  made  of  the  modelling  of  strictly  spatial  variabil- 
ity. This  is  close  to  being  accomplished  with  radar  pictures,  now  routinely  taken 
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Table  7.  Summary  of  the  Cumulative  Relative  Frequencies  and  the  Resulting  Values  of  z.  (Model  B) 
of  New  England  January  1-hr  Precipitation  in  21  Years  (1952-1972).  Areas  1 to  4 are  defined  in 
Table  3;  the  symbols  are  defined  in  the  text 
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by  the  National  Weather  Service.  Radar  scans  the  Precipitation  within  a radius 
of  125  nmi  or  250  nmi  in  a matter  of  seconds. 

A typical  picture  of  a radar-scope  survey  of  the  precipitation  surrounding  a 
station  shows  echoes  that  look  like  white  clouds  on  a black  background  (for  example, 
the  negative  of  Figure  23 >;  such  pictures  permit  an  examiner  to  judge  the  direction 
and  distance  of  storms  from  the  station.  Time  changes  give  the  direction  and 
speed  of  movement,  making  the  radar  primarily  a tool  of  short-range  forecasting. 
But  pictures  taken  regularly  at  many  stations  in  the  United  States,  for  six  years 
or  more,  provide  a basis  for  a climatological  study  of  the  areal  extent  of  precipi- 
tation. The  results  of  such  a study  are  reviewed  and  interpreted  below. 

The  outer  circle  in  Figure  23  has  a 
radius  of  125  nmi,  the  innermost  circle 
has  a radius  of  25  nmi,  and  the  second 
last  circle  has  a radius  of  100  nmi.  Be- 
tween the  25  nmi  and  100  nmi  circles, 

2 

the  area  of  29452  nmi  is  divided  into 
2 

64  cells,  each  460  nmi  . Two  adjacent 
cells  of  the  innermost  ring  have  twice 
the  single-cell  area,  four  cells  have 
four  times  the  single-cell  area,  the 
whole  inner  ring  has  eight  times  the 
single-cell,  two  rings  have  16  times 
the  area,  three  rings  have  32  times  the 
basic  area,  and,  lastly,  the  whole  us- 
able space  has  64  times  the  basic  area. 

This  made  it  possible  to  study  the  fre- 
quency with  which  fractions  of  an  area 
were  covered  by  precipitation  echoes, 
as  a function  of  areal  size.  This  was 
done  at  five  widely  scattered  stations 
within  contiguous  United  States  for  the 
four  midseason  months,  using  eight 

observations  per  day  at  3-hr  intervals.  Six  years  of  radar-picture  records  were 
used  (1969-1974).  At  each  station,  in  each  month,  therefore,  there  were  a total  of 
1440  or  1488  radar  pictures  surveyed.  On  the  average,  there  were  less  than  119 
missing  pictures  per  month  per  station. 

As  an  example.  Table  9 shows  the  cumulative  relative  frequency  of  the  frac- 

2 

tional  coverage  of  each  area  from  one  cell  (460  nmi  ) to  all  cells  or  rings 
o 

(29,  452  nmi  ).  The  kind  of  information  in  Table  9 is  ready  for  entry  in  Figures  9 
and  11  through  19.  The  frequencies  of  the  first  line,  for  10/10  coverage,  are  used 


Figure  23.  Sample  of  a PPI  Radar- 
Scope  Picture  of  Storms  Within  125  nmi 
of  the  Radar  Station.  The  area  between 
the  circles  of  radii  25  nmi  and  100  nmi 
has  been  divided  into  64  cells,  each  of 
460  nmi^ 
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Table  9.  The  Cumulative  Relative  Frequency  of  Each  Fractional  Coverage  of  Areas  of  Increasing  Size 
by  Radar  Echoes.  The  example  is  for  Cape  Hatteras,  NC  July,  based  on  6 years  (1969-1974)  of  radar 


r 


1 


I 


with  Figure  9,  those  on  the  second  line  for  9/10  coverage  are  used  with  Figure  19, 
•••,  those  on  the  line  for  1/10  coverage  are  used  with  Figure  11, 

Since  the  single-point  frequency  of  a radar  echo  was  0.  04  (Table  9),  the  cum- 
ulative relative  frequency,  0.90,  was  entered  on  the  left-hand  side  of  Figure 

9 and  the  interpolated  curve  was  followed  till  the  value  Pj.  = (1-0.0045)  was  reached. 
(The  reading  0.0045  is  on  the  right-hand  scale.)  The  abscissa  showed  a value  of 

z 4.  20.  Likewise  a z-value  was  found  for  each  pair  of  values  of  P and  Pr 

o f 

(Table  10),  ultimately  leading  to  the  geometric  mean  value  of  the  scale  distance, 
r 0.79  nmi  or  1.4(1  km.  (The  z-values  were  not  read,  or  recorded,  in  instances 
when  the  curves  in  Figures  11  through  19  were  too  flat  for  a reliable  reading.  ) 

Once  the  scale  distance  is  determined,  or  assigned  a value,  it  can  be  used 
with  Figures  9 and  11  through  19  in  a reverse  procedure  to  estimate  the  probability 
of  fractional  cover,  by  precipitation,  or  areas  of  varying  size.  Using  the  value 
r 0.  79  nmi,  this  was  done  for  Cape  Hatteras,  July  rainfall,  whose  single-point 
probability  is  0.04.  The  resulting  estimates  are  shown  (Table  11)  to  compare  with 
the  6-year  sample  frequencies  (Table  9). 

For  the  four  midseason  months  at  five  stations  the  estimates  of  scale  distance 
(r)  were  as  shown  (Table  12),  together  with  geometric  means  by  station  and  by 
month.  Table  12  shows  large  variations  in  the  estimate  of  scale  distance,  which 
make  it  difficult  to  associate  its  magnitude  with  geography  or  climatic  regime. 
There  may  be  a lesser  horizontal  persistence  at  Key  West,  Florida,  at  a latitude 
lower  than  25°N,  compared  with  that  at  other  stations,  all  of  which  are  above  35°N. 
But  this  apparent  difference  could  be  due  to  sampling  errors  and  noise  in  the  radar  - 
scope  data. 

There  is  a reasonable  variation  with  season  of  the  year.  The  horizontal  per- 
sistence or  the  scale  of  weather  patterns  should  be  greater  in  winter  than  in  sum- 
mer, and  this  is  reflected  in  the  relative  sizes  of  the  mean  scale  distance:  3.  63  km 
in  January  versus  1.  48  km  in  July. 

Making  r = 1.48  km  as  the  universal  scale  distance  in  July  at  all  seasons, 
estimates  were  made  of  the  probability  of  each  fractional  coverage  by  radar 
echoes  at  each  of  the  five  stations  in  the  sample,  to  compare  with  the  correspond- 
ing G-year  sample  frequencies.  The  overall  rms  difference  was  0.0067,  which  is 
less  than  1 percent  >'n  the  probability  of  areal  coverage.  But  the  rms  difference 
varies  systematically  from  a low  of  0.  12  percent  in  the  probability  of  full  coverage 
to  a high  of  1.  2 percent  in  the  probability  of  1/10  areal  coverage.  Figure  24  shows 
graphically  the  difference  between  the  Model  B estimates  of  the  probabilities  of 
3/10  areal  coverage  and  the  sample  frequencies  (marked  by  the  x's)  at  Key  West 
where  the  single-point  probability  is  0.02  and  at  Cape  Hatteras,  North  Carolina 
where  the  single-point  probability  is  0.04.  The  model  estimates  of  probability 
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Table  12.  The  Sample  Estimates  of  Scale  Distance  (r)  for  Each  of  Five  Stations  in 
the  Four  Midseason  Months 


■ “ 

r 

January 

April 

July 

October 

Annual 

Geometric 

Mean 

Station 

(nmi) 

(nmi ) 

(nmi) 

(nmi ) 

(nmi ) 

km 

Minneapolis,  MN 

1.42 

1.25 

1.31 

1.60 

1.39 

2.57 

Key  West,  FL 

1.54 

1.33 

0.47 

0.81 

0.94 

1.74 

Wichita,  KS 

2.72 

1.59 

0.87 

1.32 

1.49 

2.75 

Cape  Hatteras,  NC 

2.41 

1 .69 

0.79 

0.84 

1.28 

2.37 

Evansville,  IN 

2.01 

1.45 

0.76 

1.96 

1.44 

2.67 

5-Station 
Geometric 
Mean  (nmi) 

1.96 

1.45 

0.80 

1.23 

1.29 

Geometric  Mean  (km) 

3.63 

2.69 

1.48 

2.28 

2.39 

Figure  24.  An  Adaptation  of  Figure  13  to  Show  the 
Estimates  of  Probability  (solid  curves)  of  3/10 
Areal  Coverage,  by  Echo-Producing  Rainfall  in 
July  in  the  Vicinity  of  Key  West,  Florida  and  Cape 
Hatteras,  North  Carolina.  The  horizontal  scale 
has  been  converted  to  area  (km2).  The  X's  are  the 
six-year  sample  estimates  to  compare  with  the 
model  estimates 
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could  well  be  favored  over  the  noisy  sample  estimates.  Model  B implies  that  the 

Key  West  sample  of  six  Julys  underestimates  the  probability  of  3/10  areal  coverage 

2 2 

by  1 percent  or  less  for  all  areas  from  1000  km  to  50,  000  km  . 

The  model  extension  of  the  above  results  is  illustrated  further  in  Figure  25. 

If  the  single-point  frequency  is  0.07,  then  Model  B estimates  no  probability  higher 
than  25  percent  for  1/10  of  an  areal  coverage.  When  the  scale  distance  is  1.  48  km, 
this  probability  is  highest  for  an  area  of  100,  000  km".  For  greater  areas  it  will 
decrease  eventually  to  zero. 

Again,  if  the  single -point  probability  is  2 5 percent,  then  the  probability  of 
1 10  of  an  areal  coverage  will  logically  increase  continuously  with  the  size  of  the 

area.  But  Model  B specifies  this  probability  as  exceeding  50  percent  in  an  area 
2 

exceeding  5000  km  . 
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Figure  25.  An  Adaptation  of  Figure  11  to  Show  the 
Estimates  of  Probability  (solid  curves)  of  1/10 
Areal  Coverage  by  Echo -Producing  Rainfall  in  July 
When  the  Single-Point  Probability  is  0.  02,  0.  04, 

0.  07,  or  0.  25.  The  X's  are  the  six -year  sample 
estimates  at  Key  West,  Florida  and  Cape  Hatteras, 
North  Carolina 


».  MMMWB  WII  COM  I I r>IO\S 

The  stochastic  models  described  in  this  paper  and  the  resulting  applications 
(primarily  in  Figures  6 through  9,  11  through  19)  have  been  demonstrated  on  sev- 
eral practical  problems  in  the  areal  extent  of  weather  conditions  or  events.  Two 
models  were  developed:  Model  A preferred  for  macroscale  events  and  Model  B for 
mesoscale  events.  Their  limited  success  is  a strong  recommendation  for  further 
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efforts  at  improved  modelling.  There  is  the  compelling  realization  that  the  above 
results  have  been  dependent  on  Monte  Carlo  simulation.  The  challenge  remains 
to  solve  the  problem  analytically,  to  present  the  probabilities  through  algorithmic 

equations  that  could  be  used  in  computer  programming,  as  was  done  for  the  dura- 

1 2 

tion  of  n-day  precipitation  by  Todorovic  and  Woolhiser. 

Compared  with  the  direct  estimates  of  cc.  Model  A has  overestimated  the 
horizontal  persistence  of  100-mb  temperatures  over  long  distances,  while  Model  B 
has  underestimated,  somewhat,  the  horizontal  persistence  of  surface  precipitation. 
But,  for  the  intended  purpose  of  estimating  the  probability  of  a full  or  fractional 
areal  coverage,  the  models  demonstrate  a workable  simplicity,  requiring  a single 
parameter  (r)  for  each  season  but  apparently  not  for  each  geographic  regime. 

The  parameter  of  Model  B that  measures  the  horizontal  persistence  and  is 
termed  the  scale  distance  (r)  shows  a logical  decrease,  in  wintertime,  from 
9.45  km  for  24 -hr  precipitation  to  5.02  km  for  1-hr  precipitation  to  3.  63  km  for 
the  "instantaneous"  precipitation  rate  of  radar  echoes.  Likewise,  in  summertime 
the  scale  distance  has  decreased  from  5.  19  km  for  24-hr  precipitation  to  3.  61  km 
for  1-hr  precipitation  to  1.48  km  for  an  "instantaneous"  rate. 

When  the  scale  distance,  r(km),  is  known.  Figures  6 through  9,  11  through  19 
give  the  desired  probability  estimates  to  correspond  to  a given  single-point  prob- 
ability. At  a selected  value  of  z on  the  horizontal  scale  the  corresponding  area 
2 2 

s'  (km  ),  is  obtained  from 
s'  = r • 2Z 

2 

If  the  area  is  given  as  s'  , or  a line  of  travel  is  given  as  s',  then  the  z-value  for 
entry  into  Figures  6 through  9 and  11  through  19  is  given  by 

z = (in  s'  - In  r)/tn  2 


12.  Todorovic,  P.  and  Woolhiser,  D.A.  (1975)  A stochastic  model  of  N-day 
precipitation,  J.  Appl.  Meteorol. , 14:17-24. 
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